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Abstract 

This paper proposes an optimized sensing based channel access framework for the LTE cogni¬ 
tive femto-cells, with an objective of maximizing the femto-cells downlink throughput. Cognitive 
femto-cells opportunistically transmit on the macro-cell channels when they are free of use. Those 
free channels are located by means of spectrum sensing using energy detection. Moreover, periodic 
sensing is adopted to detect any changes of the sensing outcomes. The maximum attainable femto- 
cell downlink throughput varies with the macro-cell channel occupancy statistics. Therefore, the 
LTE macro-cell occupancy is empirically modeled using exponential distributions mixture. The 
LTE cognitive femto-cell downlink throughput is maximized by compromising the transmission 
efficiency, the explored spectrum opportunities and the interference from the macro-cell. An 
analytical solution for the optimal periodic sensing interval that maximizes the throughput is 
found and verified by simulations. The obtained results show that there is indeed a single periodic 
sensing interval value that maximizes the LTE cognitive femto-cell downlink throughput. At the 
peak of the macro-cell traffic, our framework increases the femto-cell throughput by ~ 15% 
compared to the senseless case. The impact of the available number of channels for opportunistic 
access is studied and no significant impact is found for more than three channels. 

Index Terms: Two tier LTE network, Cognitive femto-cell, Channel occupancy. Downlink 
throughput. Periodic sensing. Energy detection. 

L Introduction 

Mobile operators are facing an explosive growth in data traffic which imposes a great 
challenge on how to be handled with the available radio resources. Many studies have shown 
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that this huge data traffic is mostly indoor originated pQ. Therefore, a promising solution for 
this challenge is to approach towards more distributed networks architecture using low power 
short range indoor access points. In this regard, the third generation partnership project 
(3GPP) provides standards for next generation cellular systems that support deployment 
of reduced scale plug and play access points. These access points are known as femto-cell 
base stations which are connected to the mobile core network through internet cloud. [2]. 
Femto-cell base stations together with the macro-cell base station (MBS) serving the same 
geographical area form a network topology termed as two tier heterogeneous network. 

From radio network planning prospective, femto-cell base stations can share the spectrum 
owned by the MBS on an opportunistic base under cognitive radio (CR) framework [3]. 
With CR, secondary unlicensed users can access the radio channels when they are not 
being utilized by their primary licensed users. These free radio channels are called spectrum 
opportunities [1]. Therefore, letting the femto-cell base stations to take the role of secondary 
users and the MBSs to act as primary users, would facilitate achieving higher network 
throughput with more efficient use of the licensed spectrum owned by the network operator 
[S]. Femto-cell base stations sharing the spectrum with MBSs are called cognitive femto- 
cells. The potentials of cognitive femto-cells are extensively investigated in the literature, 
for example, see i, 0,0. As a 3GPP standard, LTE femto-cell base stations deployment 
is supported which will be referred to as FBS hereafter. In this paper, an optimized channel 
access framework for the FBSs that maximizes their throughput is proposed. Below are 
the connections between what has been done in the literature and the contributions of this 
paper. 

To the best of the authors’ knowledge, it has always been assumed that an accurate 
spectrum sensing is performed by the FBS [5], [S]- Based on that assumption, FBSs are 
considered to continuously have perfect knowledge regarding the spectrum opportunities 
availability. This assumption is reasonable if the concern is the instantaneous status of a 
specific channel at sensing especially with considering the high transmission power of the 
MBS that eases their detectability. However, there is no guarantee that this status will hold 
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to the end of using the channel by the FBS or to the next sensing as the MBS may resume its 
transmission. Consequently, MBS activity statistics directly influence the interfered portion 
of FBS transmission. Hence, obtaining these statistics of MBS activity is an important 
prerequisite for system’s parameters setting. In this regard, many studies have been carried 
out to characterize the downlink cellular channel occupancy statistical distribution. In pro], it 
is shown that exponential distribution can approximate mobile telephony channel occupancy 
with its advantage of traceability in hnding analytical solutions for optimization problems. 
Therefore, exponential distribution has been intensively used to model cellular channel 
occupancy, see nn, [El, na as examples. Nevertheless, many empirical studies have shown a 
heavy tail behaviour for the cellular channels occupancy which is not properly characterized 
by exponential distributions [T3|. Accordingly, some other heavy tail distributions are used 
as alternatives to model the cellular channel occupancy, among which are log-normal, jl5j . 
|16j . Beta and Kumaraswamy distributions [m. Aiming at preserving the advantage of 
exponential distributions traceability and characterizing the heavy tail behaviour of cellular 
channel occupancy at the same time, in this paper, exponential distributions mixture is used 
for LTE MBS downlink channel occupancy modeling. Using of exponential distributions 
mixture is firstly used for tele-traffic modeling in |TS| and adopted in [12] to model channel 
occupancy in WiFi systems. Moreover, using exponential distributions mixture is motivated 
by the fact that combining multiple distributions linearly improves the goodness of fitting 
compared to the case of single distribution fitting. 

Primary system channel occupancy model based optimization of sensing parameters has 
been an attractive research problem in CR arena. In HU the objective of the sensing 
parameters optimization is set as to maximize the transmission efficiency subjected to 
interference mitigation. In [T2], mutual iterative optimization algorithm for sensing time 
and periodic sensing interval is introduced. The algorithm presented in [12] optimizes the 
sensing parameters from sensing reliability prospective and no interference or throughput 
considerations are taken into account. In na, a multi-channel system is considered for 
periodic sensing intervals optimization. The driving objectives in [TB] are to maximize the 
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utilization of the available spectrum opportunities and to minimize the idle channel search 
delay. Sensing time is optimized in [20] with an objective of secondary system throughput 
maximization constrained by primary system protection. In 121 . no channel occupancy 
model is considered. Instead, the instantaneous throughput when a free channel is available 
is optimized. In contrast to di, na. la. m, this paper tackles the sensing parameters op¬ 
timization from different prospectives and addresses some of the simplihcations and general 
assumptions carried out in the previous related work. At hrst, a specihc scenario of dynamic 
spectrum sharing is considered, that is, LTE FBSs accessing the spectrum used by the MBSs 
opportunistically. By assuming this specihc scenario, the general assumptions which affects 
the optimization procedure are eliminated. Moreover, more accurate results are obtained by 
using the actual transmission and deployment parameters for the specihc sharing system. 
For example, in nn, na, na, the starting point is an assumption of a specihc value of 
the received signal-to-noise-ratio (SNR). In this paper however, the deployment scenario, 
transmission power and the propagation models for LTE systems specihed in the 3GPP 
standards are used. Furthermore, a simplihed model of exponential distribution for the 
licensed user channel occupancy is used in mi. III. 1131 while, in this paper, the statistical 
model for the channel occupancy is empirically developed, i.e., mixture of exponentials 
distributions. Furthermore, another added value for this paper is a study of the changes 
of the throughput with a real-life traffic variations. 

To summarize, the distinct contributions of this paper are 

1) Building an empirical statistical model for a primary system to facilitate the secondary 
system channel access optimization. 

2) Considering a specihc opportunistic spectrum access scenario of FBSs sharing spectrum 
with MBSs. 

3) Proposing and evaluating an FBS throughput maximization based spectrum sensing 
framework. 

The rest of this paper is organized as follows: Section [TTl presents the system model 
including deployment model, opportunistic spectrum access model and energy detection. 
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The empirical model for the downlink LTE channel occupancy is presented in Section IIIII In 
Section HVl the throughput based sensing optimization framework is explained. The obtained 
numerical results with their interpretations are included in Section Finally, Section IVII 
provides the concluding remarks of the paper. 

II. System Model 

In this section the opportunistic spectrum sharing model adopted by the FBS is intro¬ 
duced. Besides, the two state Markov channel occupancy model is explained. Moreover, 
energy detection (ED) preliminaries are presented in this section. 

A. Deployment and sharing model 

The FBS opportunistically accesses the same downlink spectrum assigned for the MBS 
which is divided into L channels. The FBS is responsible for detecting the free channels 
available for its usage with no assist from the MBS. Moreover FBS checks for MBS trans¬ 
mission reappearance using periodic sensing. In case of MBS active transmission found, FBS 
releases the channel and starts looking for another free channel. Otherwise, FBS operation 
on the same channel continues. Due to this discrete sensing, there are occasions of FBS and 
MBS simultaneous transmission. This concurrent transmission is assumed to result in an 
interference from the MBS into FBS but not in the other direction under the assumption 
that MBS have considerably higher transmission power than FBS. Moreover, the lengths 
and amount of the instances of mutual transmission on a same channel are determined by 
the instances between the sensing occasions known as periodic sensing interval. However, 
assuming that the FBS can perform one task at a time, either sensing or communicating, 
then periodic sensing is an overhead that reduces the time during which the FBS can utilize 
free channels. Accordingly, Ending the optimal periodic sensing interval is an optimization 
problem which can be objectively solved. The objective followed in this paper to handle this 
optimization problem is the FBS throughput’s maximization. FBS throughput is chosen as 
an objective as the ultimate goal of FBS deployment is to provide high data rate indoor 
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Fig. 1; Deployment and spectrum sharing model. 


cellular link. Another assumption considered in this paper is the inter-FBSs interference 
which is assumed to be small and can be considered as a part of the background noise. 
Furthermore, due to their plug and play characteristics, it is reasonable to assume that 
FBSs have high duty cycle when they are active compared to the MBS. Therefore, for 
simplihcation and due to absence of FBS traffic models, one can assume that an FBS is 
actively transmitting as long as it is powered on. Fig. [T] depicts the spectrum sharing model 
followed in this paper. 

B. Opportunistic spectrum access model 

Hereafter, the term channel will be used for an LTE MBS channel. From occupancy 
point of view, the channel can be modeled as a two state Markov process. These two 
states are: ON state representing occupied channel state and OFF state when the channel 
is idle. Each of ON and OFF states temporal length is a random variables (RV) with a 
specific statistical distribution. ON and OFF temporal length are assigned the RVs x and 
y respectively throughout the rest of this paper. The statistical distribution of x and y will 
be discussed in details with an empirical modeling in Section HTTl Channel utilization factor 
or duty cycle, u, is defined as the fraction of time during which the channel is being utilized 
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Fig. 2: Opportunistic channel access model. 


by its primary user, i.e., MBS in this paper, u is mathematically obtainable as 


E{x} 

E{x} + E{yy 


( 1 ) 


with E{y denoting the expected value. The FBS senses a channel and starts to use it in 
case of MBS transmission absence. Otherwise, the FBS senses another channel. The sensing 
is performed periodically with a period of T either to look for a free channel if none was 
found in the previous sensing or to detect MBS operation reappearance on the channel being 
used by the FBS. Transmission efficiency, r] is dehned as the ratio between the time spent 
is utilizing a free channel, T, and the time spent on both utilizing the channel and sensing 
it. Therefore, r] is found as 


T 

T -\- tg 


( 2 ) 


where tg is the time required to perform the sensing called sensing time. 

As the sensing is performed periodically in a discrete points in time, the following two 
situations are experienced 

1) The state of channel I changes from ON to OFF state one or more times within a 
period of T, meanwhile, the FBS captures a fraction of opportunities, call it C and 
misses (1 — C) of the opportunities on that channel. A formula for hnding C is provided 
in Subsection IIII-BI and Appendix A. For a system with L channels, assume that the 
captured opportunities on channel I is denoted as Q. Therefore, the whole system 
















captured opportunities, Cs, is obtained as 

c.=i-n(i-o)^ (3) 

i=i 

2) The state of channel I changes from OFF to ON state one or more times within a period 
of T while the FBS is utilizing it. Therefore, during a fraction of T both the MBS and 
the FBS use the same channel. This fraction of time of mutual operation is derived in 
Subsection IIV-AI 

The opportunistic channel access model for FBS with the captured/ missed opportunities 
and the mutual operation is depicted in Fig [2] where the ON states are represented by the 
higher level of binary representation and the OFF states are represented by the lower state. 

C. Energy detection sensing 

FBS performs spectrum sensing to locate a free channel. Spectrum sensing is basically the 
process of testing the existence of a primary user signal within a specihc band or channel. For 
spectrum sensing, there are many proposed techniques in the literature which are surveyed in 
j21j . [22] . [22] . Among the proposed spectrum sensing techniques, ED is the simplest where 
received signal energy is compared with the background noise and accordingly primary 
user signal existence or absence is declared. The price of the ED simplicity is a degraded 
performance in low SNRs. However, in cellular systems and particularly in urban areas, the 
downlink SNRs are strong due to the high transmission power of the base stations. Hence, in 
accordance with this assumption, ED can be used as a reliable technique by FBS to detect 
the free channels. 

When a channel is sensed, one of two hypotheses is declared: either the channel is free 
which is called the null hypothesis, "Ho, or the channel is used which is referred to as the 
positive hypothesis, "Hi. For the null hypothesis, noise only exists in the channel while for 
the positive hypothesis, noise bearing signal is received. Mathematically, the two hypothesis 
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framework is expressed as 


r{t) = 


n{t) 


-Ho 


(4) 


^ r(m)p< p 


m=l 


Otherwise 

H 


s{t) + n{t) 

where r{t), n{t) and s{t) are the received time domain signal, noise only components and 
signal components respectively. Accordingly, ED process is expressed as 

M 

\r(m)\‘^< p 'Hn 

(5) 

where m is an integer index for an M discrete samples for the received signal, p is called 
the detection threshold which is obtained as explained later in this section. Assume that 
implementing (jS]) is dominated by collecting the samples. Then by using Nyquist sampling 
criterion, the M samples are collected in a time of tg where M = 2Bts, with B denoting the 
channel bandwidth. 

Two statistical performance metrics to evaluate spectrum sensing techniques are used, 
namely: the probability of false alarm, pja, and the probability of detection, p^. The prob¬ 
ability of false alarm is the probability of wrongfully declaring signal existence while noise 
only is present. The probability of detection is the probability of detecting a signal that 
correctly exists. The probabilities of false alarm and detection are derived in m , jI 2 ] as 


^p-2tsB{^o + l)c^r 


Pd = u-Q 


2\/t7S(7o -l- 1) 




( 6 ) 

(7) 


where Q{-) is the Q function representing the complementary cumulative distribution func¬ 
tion (CCDF) of a Gaussian random process, cr^ is the noise variance and 70 is the received 
SNR. If the energy detector is set to achieve a specihc probability of false alarm, then the 
detection threshold, p is calculated as 


p = 2\/tgBalQ ^ 




( 8 ) 
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where Q ^(•) is the inverse Q function. 

III. Empirical LTE downlink channel Occupancy modeling 

This section presents the MBS channel occupancy statistical model developed as a part of 
this paper contribution. The section starts with introducing the measurements methodology 
and setup followed by the procedure of exponential distributions mixture htting. At the end 
of the section, the htting results are shown. 

A. Measurements 

To model the LTE downlink channel occupancy, a measurement campaign was held to 
fetch the empirical data. The measurements were held in an indoor location in Kista, 
Stockholm, Sweden in a university campus located at the GPS coordinates: 59° 24' 19.13 ” N, 
17° 56” 56.12” E. The measurements location is located in one of the industrial centres of 
Stockholm where the surroundings are densely occupied by offices, a shopping mall and 
residential buildings. Fig. [Ha] shows a Google map of the measurement location. The MBS 
signal is captured by an omnidirectional antenna connected to a real time spectrum analyzer 
(RTSA) to obtain a real time measurement data. The RTSA is controlled through a computer 
which is also used to record the data for further analysis. Fig. [3b] shows the measurements 
setup. The measured MBS signal is 40 MHz width (the full RTSA bandwidth) centred at 
2650.0 MHz. To cope with the variation of the traffic the collected measurements are divided 
into periods of four hours each. The RTSA measures the signal in the whole bandwidth 
every 30.0 ms. The 40 MHz band is composed of different LTE channels with different 
bandwidths. Hereafter, the hndings for a 5 MHz channel lies between 2653.4 and 2658.4 
MHz are considered. The ON and OFF states are found using ED with a probability of false 
alarm of 1 x 10“°. The received signal is found to have an SNR of 70 > 16 dB. Therefore, 
there is no practical risk of miss detecting some of the signals. After performing the detection, 
the ON and OFF temporal lengths are htted to exponential distributions mixture with the 
methodology explained in the following subsection. 
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Fig. 3: (a) Measurements location and (b) Measurements setup. 


B. Exponential distributions mixture fitting 

The idea of exponential distributions mixture fitting is to fit different parts of the distri¬ 
bution into different independent exponential distributions linearly combined with different 
weights. Mathematically, for a random variable, 6, the exponentials mixture distribution 
probability density function (pdf), f{0), is expressed as 

k 

f{9) = J2wAe-^^^, (9) 

i=l 

where Wi is the weight of the distribution, Aj is its corresponding exponent parameter 

k 

and = ^- 

i=l 

Following is the procedure of hnding the distribution parameters as the essence of |18j . 
The values of A* can be chosen to satisfy Ai < A 2 < ... < Afc,Vh By this, the distribution 
characterized by the exponent parameter Ai will decay the slowest compared to the other 
distribution while the one characterized by Xk will decay the fastest. Therefore, we can 
assume that some part of the tail is exclusively characterized by the hrst exponential 
distribution. Accordingly, 


1 — wie = F’(ci), 


(10a) 
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and 


k 

= 0 for 0 > Cl, ( 10 b) 

i=2 

where ci is a sufficiently large value located in the tail of the empirical cumulative distri¬ 
bution function (CDF), F{6). In the same way, if is a constant that satishes b > 1, then 



1 - = F{bci), 

( 11 a) 

and 

k 

'^Wie~^"^ = 0 ior 9 > bci. 

i=2 

(lib) 

Subsequently, Ai 

and Wi are obtained by 



^ 1 ^ 1 - F(ci) \ 

^ {b - l)ci yl “ ) ’ 

( 12 a) 


nil = (1 - F(ci))^^"^ . 

(12b) 

After hnding Ai 

and Wi the pairs (A*, Wi) ioi 2 <i < k are found recursively as 



, 1 / 1 - F(Ci} \ 

* {b-l)c,^\l-F,{bci))' 

(13a) 


Wi = 1 — Fi{ci)e^^^\ 

(13b) 


where 

Ci = a > b, 

i=i 


F,{e) = F{e). 
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Finally, the last pair {\k,Wk) is found as 


fc-i 


Wk = i- Yl 


w 


3-’ 


\k = —In 


i=i 

Wk 


(14a) 

(14b) 


Cfc yl Fki^Ck) j 

Setting the values of Ci, 6, and a is explained in [TS] . 

As stated in Subsection III-Bl ON and OFF temporal lengths are assigned the RVs x and 
y respectively. Suppose that both x and y are htted to k exponential distributions mixture 


as 




i=l 


fiy) 


(15a) 

(15b) 


i=l 


where 

wf is the weight for the exponential distribution number i for the RV x, 

Af is the exponent parameter for the exponential distribution number i for the RV x, 
tcf is the weigh for the exponential distribution number i for the RV y, 

Xf is the exponent parameter for the exponential distribution number i for the RV y, 

From the exponential distributions parameters, the channel utilization factor u is found as 


k wf" 

Si 


u = 


k ( Wf\ k (w\ 

(^j ‘‘'iS 


(16) 


With a periodic sensing performed each T seconds, the captured opportunities C are found 
using 

/A -j/i* , 

(17) 




-\^T 


Proof: see Appendix A. 
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C. Measurements and fitting results 

To investigate the advantage of using exponential distributions mixture over some other 
heavy tail distributions, the measurements hndings for x and y are htted for different 
distributions including exponential, log-normal, generalize Pareto and exponentials mixture. 
The empirical and fitted distributions for y in the period Tuesday, October , 01, 2013, 12:00 
- 16:00, are shown as samples of the results in Fig. 01 Choosing the OFF periods and the 
time is done arbitrarily as a representative case of the results. As Fig. 0] shows, htting with 
a mixture of eight exponential distributions outperforms all other fitted distributions. The 
goodness of £t is quantitatively evaluated using the log-likelihood estimation defined as [21] 



(18) 


where ^{6\f{6)) is the log-likelihood estimation for a random variable 6 having an empirical 
pdf f{6) htted to a distribution with a pdf g{0). In the case of exponential mixture htting, 
^{9\f{9)) is found as 



(19) 


Table D shows the log-likelihood estimation for diherent distributions. The log-likelihood 


estimation is calculated for htting both x and y in the same period as the one used to 
generate FigJU 


IV. Throughput driven Sensing 


A. MBS-FBS Interference Model 

Taking into account that the FBS has a very low transmission power compared to the 
MBS, the interference generated by the mutual operation is assumed to be in one direction, 
from the MBS to the FBS. Generally, in the spectrum sensing based secondary access 
framework, depending on the sensing hnding, the mutual operation can happen in three 
cases exhibited in Fig. 0]and described below. 
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Fig. 4: Empirical and fitted CDF for different distribntions. The data is for the OFF periods 
length in the period of Tnesday, October , 01, 2013, 12:00 - 16:00. 


TABLE I: Log-likelihood for htting different distribntions. The htted date is for the ON 
periods length, x, and the OFF periods length, y, in the period of Tuesday, October , 01, 
2013, 12:00 - 16:00. 


Distribution 

<F(x|/(x)) <l>(y|/(y)) 

Exponential 

0.368 0.377 

Lognormal 

0.017 0.013 

Generalized Pareto 

0.032 0.029 

k = 2 

0.292 0.244 

k = 3 

0.220 0.228 

k = A 

0.136 0.131 

Exponentials mixture k = 5 

0.068 0.066 

k = 6 

0.064 0.063 

k = 7 

0.020 0.019 

k = 8 

0.005 0.004 


1) Case 1: The secondary user miss-detects the primary signal which is active for a time 
longer than T. Accordingly, the secondary user starts to use the channel simultaneously 
with the primary user. 

2) Case 2: The secondary user miss-detects the primary signal as in case 1. However, the 
primary user changes its activity within a period of T one or more times. 
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Fig. 5: Mutual operation cases between primary and secondary users. Case 1 and 2 are due 
to miss detection while case 3 is due to primary transmission reappearance. The binary 
representation is for the primary user channel occupancy. The solid lines plot represents a 
detectable signal while the dashed lines plot corresponds to a miss-detected signal. 


3) Case 3: An absence of a primary signal is correctly detected and secondary transmission 
took place accordingly. However, after a time less than T, the primary user resumes its 
operation and possibly changes its activity for once or more. 

Due to the high transmission power of the MBS, the probability of miss-detection is 
practically zero which eliminates the occurrence of case 1 and 2. Accordingly, the fraction 
of time for mutual operation, call it r, is computed as 

r = ;^(l - P/a) (F,(T)) uT = (1 - p/„) il - ^ . (20) 

To calculate the interference due to mutual operation between the MBS and the FBS, the 
interfering power from the MBS needs to be calculated. In j^S], a model for propagation loss 
and shadow fading in LTE systems is developed. This model calculates the outdoor-to-indoor 
loss (0-1), PLm, as 


PLm = 36.71ogio(i?) + 261ogio(/c) + 0.5d + 42.7, (21) 

with R, d and fc denoting the spatial distance between the MBS and the building containing 
the terminal, the indoor distance between the wall and the terminal and the carrier frequency 
respectively. The shadow fading is modeled as a zero-mean log-normal distributed RV with 
a standard deviation aM, Therefore, the received signal from the MBS, P^, is log-normally 
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distributed as 

lOlogio (P“) : ctm), (22) 

where A/'(/i, cr) denotes a Gaussian distributed RV with a mean of /i and a standard deviation 
of (T, jjLM = — PLm^ and is the MBS transmission power. 

To calculate the signal-to-interference ratio received at the terminal, the received power 
from the FBS needs to be obtained. Similarly to the model for the prorogation loss for the 

MBS signal, in (23], a model for the indoor path loss for the FBS signal, PLp, is found as 

PLf = 43.31ogio(d) + 201ogio(/c) + 11.5, (23) 

with a zero mean log-normal shadow fading having a standard deviation ap- Consequently, 
the FBS received power at the terminal Pf is distributed as 

lOlogio {Pf) : AA(/iF, ap), (24) 

where jj^p = (^Pf — PLp^ if P^ is the FBS transmission power. 

B. Periodic sensing optimization 

When the FBS is transmitting in the interference free periods, the achieved FBS through¬ 
put, call it Go is limited by the SNR, 70, as 

Cq = B ■ log2(l -|- 7o). (25) 

In the same way, the FBS throughput in the interfered periods, G, is 

G = P ■ log2(l 7), (26) 

where 7 is the signal to interference plus noise ratio (SINR). For the FBS throughput during 
the whole operation time, call it Caii, Go and G are weighted by (1 — r) and r respectively 
and summed up. Moreover, it should be noted that the FBS is efficiently transmitting with 
a factor of rj during a fractional time of captured opportunities, (s, (see ([2]) and ([3])). Hence 
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Call is expressible as 


Call = {vCs) {t-C+{1-t)-Co). 


( 27 ) 


Let us define x as the ratio between the drop in the average FBS throughput due to the 
opportunistic access and the interference free average FBS throughput. Accordingly, y is 
found as 



The mutual operation fractional time, r, and the throughput drop, y, are related as 


r = 1 


vCs 


— CTn 


(29) 


Proof: See Appendix B. 

As (1201) and fl2^ are equivalent, equalizing their right hand sides yields 

X = 1 - h • C (1 - a • r) (30) 


where 


/iF 


Considering that r] and Cs are T dependant as in ([2]), (ITT)) and (jS]). Then, (150]) can be solved 
for T by specifying the targeted value of the throughput drop, y. However, specifying the 
targeted throughput drop, y, provides a suboptimal solution for the periodic sensing interval, 
T, as there may be lower achievable throughput drop with another value of T. Moreover, 
there is no guarantee that there is a valid solution for all values of y as there is a minimum 
throughput drop that can be achieved depending on the MBS traffic. The optimal solution 
for the periodic sensing interval, Topt, is achieved by minimizing the throughput drop, y, as 


Topt = arg min {1 - ?] • C* (1 - a • r)} . 

T 


(31) 
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Hereafter, the highest achievable throughput corresponding to Tgpt will be denoted as Copt- 

V. Numerical results 

As stated in Subsection IIII-Al the occupancy in a span of 40 MHz in the 2.6 GHz LTE 
band is measured. Two 5 MHz channels were found in that band when the measurements 
were conducted. The two 5 MHz channels show basically the same occupancy behaviour with 
highly correlated handled traffic during the measurements time. The empirical data and the 
modeled channel occupancy using exponential distributions mixture of the two channels 
are used as an input for a simulation study. The simulation part treats the throughput 
based opportunistic sharing framework developed in this paper. Below Table UTI shows the 
simulation and model parameters ge. 

TABLE H: Simulation and model parameters 


Parameter 

Value/Description 

Pfa 

1 X 10-3 

ts 

20 ms 

jdM 

40.0 dBm 

Pf 

23.85 dBm 

Pm 

7 dB [26] 

(jp 

4 dB [26] 


— 170 dBm/Hz 

L 

2 channels 

k 

8 


For benchmarking purposes, senseless throughput, Csl-, is defined as the FBS throughput 
when no sensing takes place. Instead, the FBS randomly chooses one of the available channels 
for its transmission. The senseless approach is adopted in [26]. As a consequence for the 
senseless strategy, whenever there is an MBS transmission, there exists mutual operation 

^The value of Pf is adjusted to give an average interference free throughput of 100 Mbps. 

^The inter FBS interference is considered as a part of the background noise. 
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and interfered FBS transmission. Fig. Ea] shows the senseless thronghpnt and the optimal 
achievable thronghpnt when FBS periodic sensing interval is optimized according to fl3ip . 
The same hgure shows the available interference free opportnnities for FBS, (1 — m). More¬ 
over, the optimal throughput, Copt, is shown in Fig. [Ha] with the corresponding optimized 
periodic sensing intervals in Fig. Ihbl All the results shown in Fig. [H]are obtained using the 
measurements data taken in different periods of the day October, 02, 2013. This day has 
been randomly picked. As depicted by Fig. [HaJ the more the interference free opportunities, 
the higher the senseless throughput, the lower the optimal interference free throughput 
drop and the lower the optimal periodic sensing interval. In accordance to that, the gain in 
the achieved throughput when periodic sensing interval is optimized reaches its peak with 
the lowest available opportunities. This result reflects the necessity of applying the periodic 
sensing intervals optimization as increasing the throughput is needed more with the increase 
off the traffic which is expected to happen as the same time for both MBS and FBSs. 


o) 



Fig. 6: (a) The senseless throughput, the optimal achieved throughput, the available 
opportunities for FBS and the minimum achieved interference-free throughput drop, (b) 
Optimized periodic sensing intervals. 


Fig. [7| shows an example of the achievable throughput with different values of the periodic 
sensing interval. The shown results are for the measurements performed during the period 
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October, 02, 2013, 08:00 - 12:00. As the figure shows, the highest obtainable throughput, 
Copt, of about 85 Mbps is attained when the optimal periodic sensing interval is used. For 
the other values of the periodic sensing interval, lower or higher than the optimal periodic 
sensing interval, the throughput is lower than Copt- The figure illustrates the credibility of 
using the optimal periodic sensing interval. 
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Fig. 7: The change of the achieved throughput with the change of the periodic sensing 
interval, T. The results are for the period October, 02, 2013, 08:00 - 12:00. 


Fig. E] shows the transmission efficiency, r], and the captured opportunities, for a single 
channel. The figure illustrates the contradictory trends for the transmission efficiency and 
the captured opportunities, i.e. longer periodic sensing interval results in higher transmis¬ 
sion efficiency but more missed opportunities. The total efficiently used opportunities for 
transmission are then decided by the product of the two quantities which is shown in the 
figure. For the whole system consideration, the captured opportunities in all channels, Cs is 
then taken into account instead of the captured opportunities in each channel separately. 
Consequently, the total efficiently used opportunities determined by the product of 7] and 
(s is higher for our two channels system as the figure exhibits. 

Fig. [Hal shows the CDF of the throughput when the periodic sensing is optimized. The 
CDF curves are bounded between two limits. The upper limit is the case of MBS traffic 
absence and noise only existence, i.e. Cq. The lower limit is the senseless throughput, Csl- 
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Fig. 8: Obtained values of the transmission efficiency, r], the captured opportunities, ( and 
the utilized opportunities for a single channel and the whole system. 


The lower limit is traffic dependant. Therefore, the period of the highest traffic is taken 
as the absolute lower limit for the whole day of measurements. Apart from Co and Csl, 
the throughput CDF curves in Fig. [9a] can be divided into two regions. The first region 
is the interfering transmission plus no opportunities region. In this region the throughput 
is degraded due to one of two reasons. The first reason is a simultaneous FBS and MBS 
transmission caused by MBS transmission reappearance while FBS is being transmitting. 
The second reason is that the sensing has resulted in no opportunities to utilize by the 
FBS. The second region is the non-interfered transmission region where the opportunities 
are being used by the FBS without MBS transmission reappearance. Fig. [9bl demonstrates 
the two regions of the throughput CDF for the period October, 02, 2013, 08:00 - 12:00 as 
an example. 

The influence of the number of the available channels, L, on the achieved optimal through¬ 
put is studied and the results are presented in FiglTUlfor the period October, 02, 2013, 08:00 - 
12:00. As there existed two 5 MHz channels in the measured band, the data for higher number 
of channels is a synthetic data generated by bootstrapping the data for these two channels 
and other channels with different bandwidths. The mean and the 5,10,50,90 and 95 per¬ 
centiles of the throughput when the number of the available channels changes from 1 to 5 
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(a) 
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Fig. 9: (a) The CDF for the throughput in the cases of MBS traffic absence, sense less 
scenario and periodic sensing optimization, (b) Example of the regions of the throughput 


CDF 


are shown in Fig [101 The results when a single channel is available basically represent the 
senseless scenario which is the worst case scenario. The diversity gain in the throughput is 
then judged by the increase of the throughput with the increase of the available number 
of channels. In this regard, it is seen from Fig. [10] that increasing the number of channels 
would increase the achieved optimal throughput. However, for more than three channels the 
optimal throughput saturates. This trend is explained as follows, with the increase of the 
number of channels, the captured opportunities found as in ([3]) would increase. Yet, this 
increase in the captured opportunities impact grows less significantly with increasing the 
number of available channels till it saturates at L = 4 channels. 

VI. Conclusions 

Cognitive LTE femto-cells whose access the LTE macro-cell spectrum on an opportunistic 
basis are considered in this paper. The paper contributes in introducing a framework for 
downlink channels access for LTE femto-cells. The channel access framework is based on 
employing periodic sensing using energy detection for finding the free channels. The frame¬ 
work objective is to maximize the LTE femto-cells downlink throughput. The throughput is 
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Fig. 10: 5,10, 50, 90 and 95 percentile of the throughput when different number of channels 
are available. 


maximized by optimizing the sensing intervals in a multi-channel spectrum sharing system 
by compromising the transmission efficiency, the captured opportunities and the interference 
from the macro-cell base station. The statistics of the macro-cell channel occupancy impacts 
the optimization procedure. This channel occupancy is empirically modeled as a part of 
the paper contribution. The analytical derivations have shown that there exist an optimal 
periodic sensing interval which maximizes the femto-cell downlink throughput. Hence, a 
more efficient use of the spectrum assigned for LTE downlink transmission is achievable 
following the framework presented in this paper. 


Appendix A 
Proof of equation I fTTH 

Let Tz to be the average opportunities during a period of T. z is either x or y denoting 
if the channel is occupied or free. The sensing point can be an end or a start of an idle 
period, so in that case we use T^(t) instead of T^. Therefore, we have four possible cases, as 
illustrated in Fig. lA.ll 
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Fig. A.l: Possible channel statuses at sensing 


Let y to denote the remaining time of an OFF period when sensing is performed. The 
distribution of y is given by [27j 

l-Fy{y) 


Fviy ) = 


E{y] 


(A.l) 


In the same way, if the remaining time of an ON period when sensing is performed is denoted 
as X , then x is distributed as 

l-Fx{x) 


Fx{x) = 


E{X] 


(A.2) 


Renewal theory is applied to hnd Ty{t),Tx{t),Ty{t) and T^it). Following are the derived 
formulas. 

OO - X f/ , » 


T.(t) = 


1 - Ex[x) 
E{x] 


Ty{t — x)dx. 


(A.4) 


OO c 

fy{t) = t J fY{y)dx + J fviy) {y + f^{t - y)) dy. 


(A.5) 


Tx{t) = / fx{x)Ty{t - x)dx. 


(A.6) 


If Laplace transform is performed for flA.Sp . flA.4D . (lA.Sp . and flA.bp . then following formulas 
are obtained 

E{v} ■ T,(S) = + (1 _ Fy(s)) Us). (A.7) 
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i?{x}-T,(s) = (l-F,(s))T,(s). 


Us) = fx{s)fy{s). 


Tyis) 


/y(0) ~ /y(g) 

c2 


+ fY(s)Tx{s). 


Therefore, 

1 - Fx(s) ( friO) - fY(s) \ 

’ E{x]-s^ [i-JyWxM)' 


(A.8) 
(A.9) 
(A. 10) 

(A.11) 

(A.12) 


If we define the missed opportunities, T, as the fraction of missed OFF periods, then, T is 
expressed as T = T^, + Ta,,. representing the opportunities missed from the start of an 
OFF period till the sensing is performed which is formulated as 

= (A.13) 

On the other hand, Ta, is the opportunities missed due to changing the status of the channel 
from ON to OFF within a period of T. is found as 


T = 


U 



(A.14) 


Sine the total sum of the missed opportunities, T and the captured opportunities, Ci is 
equivalent to the total available opportunities, then 


( = (1-m)-T. (A.15) 

Therefore, with the exponential mixture of the ON and OFF periods shown in fllSap and 
(HlbD, c is found using 




(A.16) 
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Appendix B 
Proof of equation (1^ 


At the terminal served by the FBS, having the received power from the FBS, the interfering 
power from the MBS and the noise variance results in a received SINK, 7 , calculated as 

PF 


7 = 


pM ^2 ■ 
r ' n 


(B.l) 


Cellular systems are interference dominating systems where the noise variance is neglectable 
compared to the interference [2S]. Therefore, with the log-normal distribution of both 
and PF^ as in (12^ and (Ell), 7 is approximated as a log-normal distributed RV as 


lOlogio (7) : 

where = fip — hm and cr^ = -|- alj. 

In the same way, 70 is expressed as 

pF 

r 

7 o = 

Accordingly, 70 is a log-normally distributed RV as 


(B.2) 


(B.3) 


lOlogio (70) : F'i^F - ctJ.CTf). 


(B.4) 


Noting that all the quantities in flB.4D are taken in their logarithmic (dB) scale. 
From (jUD, (EH]), ([27]) and (ESj) . 

AjrRloga (1 -h 7) + (1 - 'r)Rlog2 (1 -h 70) } 

^ ' A{Rlog 2 (l + 7o)} 

^B- e{t ■ lOlogio (1 + 7) + (1 - r) • lOlogio (1 + 7 o) } 

B- AjlOlogio (1 + 7 o) } 
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In the serving area of the FBS, the valid assnmption of 7, 70 3 > 1 leads to 


1 - X = h • G 


e[t ■ lOlogio (7) + (1 - • lOlogio (70) } 

^{lOlogio (70) } 


(B. 6 ) 



T ■ Pm) + (1 - r) • (/xf - 

(hF - (yl) 


Therefore, 


r = 


I'C. 



(B.7) 
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